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Abstract 



The present work deals with the two-dimensional incompressible, laminar, steady-state boundary layer equations. 
First, we determine a family of velocity distributions outside the boundary layer such that these problems may have 
vo , similarity solutions. Then, we examen in detail new exact solutions, called Pseudo -similarity, where the external ve- 

locity varies inversely-linear with the distance along the surface (u e (x) = UacX" 1 ). The analysis shows that solutions 
exist only for a lateral suction. For specified conditions, we establish the existence of an infinite number of solutions, 
^~>, including monotonic solutions and solutions which oscillate an infinite number of times and tend to a certain limit. The 

properties of solutions depend on the suction parameter. Furthermore, making use of the fourth-order Runge-Kutta 



scheme together with the shooting method, numerical solutions are obtained. 
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1. Introduction 



£> ' In this paper we are concerned with the classical two-dimensional laminar incompressible bound- 

ary layer flow past a wedge or a flat plate [28] . For the first approximation, the model is described by 
q^ I the Prandlt equations or the boundary layer equations 

q. 

r — [', (1-1) ud x u + vdyU = u e d x u e + udyyU, d x u- t rd y v = 0, 

o. 

where (x, y) denote the usual orthogonal Cartesian coordinates parallel and normal to the boundary 
y = (the wall), u and v are the corresponding velocity components, and the constant v > is 
the kinematic-viscosity. The function u e is a given exterior streaming velocity flow which is assumed 
throughout the paper to be nonnegative function of the single variable x; u e = u e (x), and is such that 
H I u(x,y) tends to u e (x) as y — > 00. Equations (II. ID can be written in the form 

(1.2) dy^d^ - d X ip8yyip = U e 8 X U e + VdyyyTp, 

where ip is the well-known stream function defined by u = d y ip, v = —d x ip. 
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This problem with appropriate external velocity flow has been the main focus of studies of partic- 
ular exact solutions. Research on this subject has a long history, which dates to the pioneering works 
by Blasius [B] and Falkner and Skan [J3] in which the external velocity is given by 

(1.3) u e {x) = Uoo x m ( Uoo > 0). 
Their investigations lead to solutions to (1.2) in the form 

(1.4) i>(x,y)=ax a f(byx-P), a,b>0, 

where 

m+ 1 1 — m 

(1.5) a = — - — and /3 = — - — . 

The equation for / is the well-known equation obtained by Falkner-Skan [T3] 

(1.6) r + ! ^// // + m(l-/ /2 )=0 on (0,oo), 

or, if 7Ti > —1, 

(1.7) f"' + ff" + a(l-f' 2 ) = on (0,oo), 

where 

2771 

a 



777+1 

Such equations occur at wedge flows [28} p. 170]. These equations have received considerable attention 
in the literature. We refer the reader to the works of Rosenhead [27], Schlichting and Gersten [28] . 
Weyl [30] and Coppel [11] and the references therein. Note that from (1.4) it is easily verified that, 

This means that a solution ip(x,y) for y fixed is similar to the solution ip{xQ,y) at a certain xq. This 
solution is called invariant or similarity solution and the function / is called the shape function or the 
dimensionless stream function. 

The broad goals of this paper is to study equation (1.7) when taking the limit a — ► — oo. This limit 
case, corresponding to m = — 1, is that of flow in a two-dimensional divergent channel (or sink flow) 
[28; p. 170]. We prove that transformation (1.4) is much too restrictive; that is problem (1.1) has no 
similarity solution and we shall see, by rigouros arguments, that the term -flog(x) must be added to 
the expression (1.4) and the surface must be permeable with suction to obtain new exact solutions 
which are not similarity. In addition, we shall prove that there is an infinite number of solutions. 
These results are given in Section 3. Before this analysis, we shall identify, in Section 2, external 
flows, such that problem (1.2) may have similarity solutions. The main result of this section indicates 
that problem (1.2) has solutions under the form (1.4) if the external flow is of the power-law type 
(1.3). 



2. Similarity solutions 

In this section we shall obtain external flows such that the partial differential equation (1.2) 
accompanied by the boundary condition 

(2.1) lim d y tp(x,y) = d y tp(x,(X)) = u e (x), 

has a solution under the form (1.4), where a + (3 = 1. 

The main problems, arising in the study of similarity solutions, are related to the existence of 
the exponents a and (3 and to the rigorous study of the ordinary differential equation satisfied by 
the profile / which is, in general, nonlinear. For the layer equation (1.2), the classical approach for 
identifying a and (3 is the scaling and transformation group [I]. The essential idea is to seek a and b 
such that if ip satisfies (1.2) the new function ip K (x,y) = K a i(){n b x,K,y) is also a solution. 
Let ip be a stream-function to (1.2) defined by (1.4) where a + (3 = 1. Assume first that (3 / 0. We 
choose a = — %, b = 4, and define ip K (x, y) = K a ip (K b x, rey) . Hence a + b = 1, ip = ip K and 

L(iP K )(x,y) = K a+3 L(i;)( K b x,Ky) 
for any n > 0, where L is the operator defined by 

L(l/}) = dy^d^ - d X ^)d 2 y y i) ~ V^ l/) . 

According to equation (1.2) we deduce 

h{x) = K a+3 h( K b x), 
where h{x) = u e (x)d x u e (x). In particular, for fixed xq > 

h{ K b x )=K- (a+3) h(x ). 
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h(x) = x b x b h(xo 



Setting x = n b xo we infer 

Solving the equation 

(2.2) u e —^ = x~ s ~b~x^ b ~h(xo) 

yields us 

(2.3) u\(x) = Cl x 2m + c 2 , 

for all x > 0, where m = a — [3 and c\ and c 2 are constants, for (3 7^ 0, since — ^jp + 1 = 2(a — (3). 
For (3 = 0, hence a = m = 1, the new function 

i(j K {x,y) = K a ip(K,~ a x,y), 

for any fixed a 7^ 0, is equivalent to ip and satisfies 

L(i> R )(x,y) = K a L(tP)( K ~ a x,y) 

for any k > 0. Arguing as in the case (3 7^ one arrives at (2.3) with m = 1. Next, because 
limy^oo x 2m f'(yx~P) 2 = c\x 2m + C2, the function /' has a finite limit at infinity, which is unique 
and is given by c\ + C2X~ 2m . This is acceptable only for c 2 = 0. 



The above result indicates, in particular, that for a prescribed external velocity satisfying (2.3) 
the real numbers a and (5 are given by (1.5) and condition (2.3) is necessary and sufficient to obtain 
similarity solution under the form (1.4) where a + (5 = 1. However, for a general external velocity, it is 
possible to obtain an exact solution which is not similarity solution as it is seen in [21]. In this paper 
the authors considered 

u e (x) = cix 1/3 + c 2 x" 1/3 . 

A stream-function ip is given by 

(2.4) ip(x,y) = x 2 ' 3 f(yx~ 1 ' 3 ) + cx~^y, c = const, 

where / is a solution of 

(2.5) vf + aff" -(a- (3)f 2 = £, £ = const. 

3. The pseudo-similarity solutions 

In the present section we restrict the attention to the case m = — 1 and get new solutions to the 
problem 

(3.1) Qy^l y ^ - d x ^d 2 yy i, = vd 3 yyy ^ - ulx- 3 , 

subject to the boundary conditions 

(3.2) d y il>(x,Q) =u w x~ 1 , d x ip(x,0) = -v" 1 , d y ip(x,oo) =u OQ x~ 1 , 

where v w is a real number ( v w > for suction and v w < for injection), u w and u M are nonnegative 
and satisfy u w < u^ . The velocity distribution u e = ^ is found in the case of divergent channel (or 
sink flow) [281 P- 170]. The analysis of this section is motivated by the work by Magyari, Pop and 
Keller [2H] concerning a boundary-layer flow, over a permeable continuous plane surface, where the 
x-component velocity tends to zero for y large (uoo = 0). In [26] the authors showed that if m = — 1 
problem (1.2), (2.1) has no solution in the usual form (1.4). For Uoq ^ and according to Section 2, 
the function ip can be written as 



i){x,y) = y/uu^f J yx 



^oo _1 



/' 



Since v w = ^^- s/vu^j "(0) (see the Appendix) we deduce v w = for m = — 1 and the following 
oridinary differential equation for / 

r r + / /2 -i = o, 

(3.3) 

( /'(0) = C,/'(oo) = l, 

where C = — is in the interval [0, 1) and /(0) can be any real number. Since this problem does not 
contain / it is convenient to study the equation satisfied by 9 = /'; that is 

(3.4) 




The stability of equilibrium point (1,0) of (3.4) cannot be determined from the linearization. To 
analyze the behavior of the nonlinear equation (3.4) i, we observe that 



E'(t) = 0, 



where E is the Liapunov function given by 



E(t) = 
Then, for some constant c, the following 



V(i) 2 + ±0(£) 3 -#(£). 



1 

3' 



±V2U + 



1/2 



holds. The analysis of the algebraic equation of the phase path in the phase plane reveals that the 

equilibrium point (1,0) is a center. Hence, Problem (3.4) has no solution for any £ > — 1 except the 

trivial one 9 = 1 (see Fig. 3.1). 

Note that if we impose the condition 9(oo) = — 1 instead of 9(oo) = 1, which is also of physical interest, 

it is easy to see that, for any £ < 2, there exists a unique solution up to translation. This solution 

satisfies 



and we find that 



-9>(tf + \9{tf-9(t) = 2 -, 



9{t) = 2-3tanh 2 ±t/v / 2 + arctanh |(2 - C)/3) 1/2 } 




Fig. 3.1 Classification ol solutions ol 8" + 8 2 — 1 = according to #(0) and 8'(0). 
To obtain exact solutions to (3.1), (3. 2), we look for "'pseudo-similarity " solutions under the form 
(3.5) ip(x, y) = aF(x, byx~ l ), 



where a = y/vuoo and b = y^^T- Assuming F(x,t) = f(t) + H(x) one sees that (cf. Appendix) 

r r+if' + f 2 -1 = 0, 

(3.6) 

( /'(0)=C, f(oo) = l, 

and 

i?(x) = 7logx + c, 

where j = v w (vuoo) ' and c is constant. Without loss of generality we may take c = 0, since 
ijj(x,y) = s/vu^,f(t) + u w log(x) + ^JvUoqC satisfies (3.1) — (3.2) for any real number c. Note that 7 
plays the role of suction/injection parameter. 
To study (3.6) it is more convenience to consider the second ordinary differential equation 

( 9" + 7 0' + 2 - 1 = 0, 
(3.7) 

[ 0(0) = C, 0(oo) = 1, 

where < £ < 1 and 7 7^ 0. In fact, the real number 7 will be taken in (0, 00). The existence of 
solutions to (3.7) will be proved by means of shooting method. Hence, the boundary condition at 
infinity is replaced by the condition 0'(O) = d, where d is a real number. For any d the new initial- 
value problem has a unique local solution 0^ defined in the maximal interval of the existence, say 
(0, Td),Td < 00. We shall see that for an appropriate d the solution 0^ is global and satisfies 

(3.8) d (oo) = 1. 

A simple analysis in the phase plane shows that problem (3.6) may have solutions only for 7 > 0. In 
fact, the ordinary differential equation in (3.6) is considered as a nonlinear autonomous system in R 2 , 
with the unknown (9,9'), mainly 

r 9' = <p, 

(3.9) <^ 

k ^ = -7c^+l-0 2 , 

subject to the boundary condition 

(3.10) 0(0) = C, <p(0) = d. 

The linear part of the above system at (1,0) is the matrix 

V "2 -7 
The eigenvalues of J are 

Ai 

if 7 > \/8 and for | -y | < \^8, 



-7- 


-vV 


-8 




2 




-7- 


iVS" 


7 2 



-7+ Vr 



-7 + V8-7 2 



Therefore, the hyperbolic equilibrium point (1, 0) is asymptotically stable if 7 is positive and unstable 
for negative 7. In particular problem (3.5) has no nontrivial solutions if 7 < 0. If 7 > we deduce 
from the above that there exists 5 > such that for any d and Q satisfying d 2 + (£ — l) 2 < 5 2 the 
local solution 9 d is global and satisfies (3.8). In the following we construct solution to (3.6) where 
the condition d 2 + (£ — l) 2 < 5 2 is not necessarily required. For a mathematical consideration the 
parameter £ will be taken in (— 1, Vo]. The following result deals with nonnegative values of £. Let us 
consider a real number d such that 

(3.11) d 2 <2((l-^), 

where < ( < \/3. We shall see that any local solution of (3. 9), (3. 10) is global and satisfies (3.8). To 
this end we consider again the Liapunov function E(8(t), <p(t)) = 2 ( / 9 (^) 2 + 3^ 3 — $■ Along an orbit we 
have 

±E(0(t)Mt)) = -'r<p(t) 2 <o. 

Hence 

E(9 d (t),9' d (t))<E(C,d), 

for any t < T d . On the other hand, from (3.9) and (3.10), there exists to > 0, small, such that 9 d is 
positive on (0, to)- Assume that 9 d vanishes at some t\ > to and suppose that 9 d {t\) 7^ 0. Because 

E((,d)>E(e d (t),d d (t))>±e' d (t 1 ) 2 , 

for all < t < t x . we deduce \d 2 > C(l - |C 2 ), which contradicts (3.11). Therefore 9'Jtx) = 0. In 
this case we deduce from the equation of 9 d that 9 d (ti) = 1 and then 9 d is nonnegative on a some 
neighbourhood of t\. Consequently the local solution is nonnegative as long as there exists. To show 
that 9 d is global we note that 

E((,d)>±9 d (t) 2 + ±9 d (t)-9 d (t)>~, 

for all t < T d , since 9 d is nonnegative. Hence 9 d and (then) 9' d are bounded. Consequently 9 d is global. 
It remains to show that 9 goes to unity at infinity. To this end we use the Bendixon Criterion. Let T 
be the trajectory of (9 d , 9' d ) in the phase plane (0, 00) x R for t > and let T + (T) be its w-limit set at 
infinity. From the boudedness of T it follows that T + (T) is a nonempty connected and compact subset 
of (0,oo) x K ( see, for example [U p 226]). Moreover (—1,0) T + (T), since 9 d is nonnegative. Note 
that if r + (T) contains the equilibrium point (1, 0) then r + (7~) = {(1, 0)} , since (1, 0) is asymptotically 
stable. Assume that (1,0) T + (T), Applying the Poincare- Bendixon Theorem [HJ p 44] we deduce 
that r + (T) is a cycle, surrounding (1,0). To finish, we shall prove the nonexistence of such a cycle. 
We define P(9,ip) = ip,Q(9,(p) = —7^ + 1 — 9 2 ,ip = 9 d and 9 = 9 d . The function (0,ip) satisfies the 
system 9' = P(9, 93), if' = Q(9, <p). Let D be the bounded domain of the (9, </?)-plane with boundary 
r + . As P and Q are regular we deduce, via the Green-Riemann Theorem, 

(3.12) / / (d v Q + d e P) d^pd9 = I (Qd9 - Pd^) = 0, 

J Jd Jr+ 

thanks to the system satisfied by (9, (p). But d v Q + dgP = 7 which is positive. We get a contradiction. 
To complete our analysis, we shall determine a basin of the critical point (1,0). Let 

V=l(C,d) GR 2 :C>-l,^ 2 + cQc 2 -l) <l 



Let us consider a one-parameter of family of curves defined by 



E (e,ip) = ^ 2 + h 3 



e = c, 



where C is a real parameter. Note that, in the phase plane, this family is solution curves of system 
(3.4). The curve ip 2 = 29 — |# 3 + |, corresponding to C = |, goes through the point (2,0) and has 
the saddle (—1,0) (7 = 0) as its a and u>-limit sets. We note this solution curve by 7i, which is, in 
fact, an homoclinic orbit and define a separatrix cycle for (3.4). We shall see that the bounded open 
domain with the boundary Ti is an attractor set for (1,0) of system (3.9) where 7 > 0. This domain 
is given by E(9,ip) = C,9 > -1, for all — | < C < §, which is V '. As j^E < any solution, with 
initial data in V cannot leave V . By LaSalle invariance principle we deduce that for any (£, d) in V 
the u>-limit set, r + (£, d) is a nonempty, connected subset of V n {tp = 0} , (see [U p. 234]). However, 
if 9 7^ 1, <p = is a transversal of the phase-flow, so the u;-limit set is {(1, 0)} . This means that V is 
a basin of the critical point (1, 0) of (3.9). 




Fig. 3.2 A basin of attraction V of the critical point (1,0). 



4. Numerical results 



In this section numerical solutions of the boundary-value problem (3.7) are obtained by using the 
fourth-order Runge-Kutta scheme with the shooting method. We plot the dimensionless velocity 9 in 
term of the similarity variable t, for various value of the shooting parameter d. 
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Similarity variable t 



Fig. 4.1 Velocity profiles in terms of d = 0'(0) for fixed £ = 0.2 and 7 = 0.5 
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Fig. 4.2 Velocity profiles in terms of d = 8'(0) for fixed £ = 0.2 and 7 = 2.83 



5. Conclusion 



In this work the laminar two-dimensional steady incompressible, boundary layer flow past a 
moving plane is considered. It has been shown that the problem has solutions having a similarity form 
if the velocity distribution outside the boundary layer is proportional to x m , for some real number 
m. In the second part of this paper, we are interested in question of existence of solutions in the case 
where the external velocity is the inverse-linear function; m = — 1. This situation occurs in the case 
of sink flow. To obtain exact solutions the stream function ip is written under the form 



(5.1) 



i/j(x, y) = y/uuoofit) + v w log(x). 



It is shown that the ordinary differential equation satisfied by / has multiple solutions for any v v 
positive and no solution can exist if v w < 0. A sufficient condition for the existence is derived: 



(5.2) 



o-i, ^/"(°) 2 + c(^ 



1 < 



We have obtained two family of solutions according to 7 = v w (£titoo) _ ■ If 7 > V°, /' is monotonic 
and goes to 1 at infinity, but if < 7 < y/8, we have a stable spiral. The function /' oscillates an 
infinite number of times and goes to 1. So if we are interested in solutions to (3.7) such that 



we must take u v 



-K /' < 1 

and Uoo > satisfying — u^ < u w < Uqo and v w > 



wu r 



a/2 



Condition (4.1) indicates also that for the same positive value of the suction parameter the permeable 
wall stretching with velocity u w x~ 1 ,u w > has multiple boundary-layer flows. Every flow is uniquely 
determined by the dimensionless skin friction /"(0) which can be any real number in the interval 

( — J I + 2( (1 — C 2 /3), a/ 3 + 2£ (1 — C 2 /3) ) , where ( = u w u^. The case Mqo = was considered by 

Magyari, Pop and Keller [26]. The authors showed, by numerical solutions, that the boundary layer 
flow exists only for a large suction parameter (7 > 1.079131). 

The existence of exact solutions of the Falkner-Skan equation under the present condition was dis- 
cussed by Rosenhead \27\ pp. 244-246] who mentioned that these results may be obtained by rigorous 
arguments which, in fact, motivated the present work. We note, in passing, that it is possible to obtain 

solutions if the the skin friction satisfies |/"(0)| > J± + 2( (1 - C 2 /3). 



Appendix 

Let us now derive problem (3.5). We assume that the external velocity is given by u e (x) = UocX™ 1 , 
where m is not necessary equal to -1. We recall that the stream function satisfies the following equation 

(5.3) d y ^d 2 xy ^ - a^V = v&mfi + mu loX 2m -\ 
with the boundary conditions 

(5.4) d y ip(x,0) =u w x m , d x ip(x,0) = —v w x 2 ; d y ip(x,oo) =u OQ x m . 

To obtain exact solutions to (5. 3), (5. 4), we look for "pseudo-similarity " solutions under the form 

(5.5) i/j(x,y) = ax a F(x,byx~ f3 ). 



where a 



(5.6) 



m+l 
2 ' 



ft 



m— 1 



j 1 , a = ^/vu^ and b = ^J^f-- Inserting (4.4) into (4.2), (4.3) leads to 
' F'" + ^FF" - m(F' 2 -l)+x (F"d x F - F'd x F') = 0, 

F'(x,0)=C, F'(x,oo) = l, 
1 + m 



-F(x,0)+xd x F(x,0) 



^/vu c 
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where the primes denote partial differential with respect to t = x/^jf-yx 2 • By writing 

F(x,t) = f(t)+H(x), 
we find 

' /'" + ^Fff" - m if 2 - 1) + /" ( xH ' + ^h) = 0, 



(5.7) 



— s— /(°) + -^— H ( x ) + ^^ (») = -7== 



/ / (o) = Ce[o,i), /'(oo) = i. 



Hence, there exists a real number 7 such that 



(5.8) 



r + H^/r -m(V 2 -i)+ 7 r = o, *>o, 



xH' + ±±™# = 7, x > 

1 + 771 



-/(0) + 7 



2 V^- 

First, let us assume that m^-1. The solution i7 is given by 

TT/ s _i+m 27 

H[x) = ex 2 + 



, /'(0) = C, /'M = l. 



1 + m' 



where c is a constant, and then ip(x, y) = ac + ax 2 / /(£) + — 2- J . The new function g = f + jt^ 
satisfies the Falkner-Skan equation. Thereafter, we will assume that m = — 1 and this leads to 



(5.9) 



and 



en , -f/2 



r+if+f -1 = 0, 

/'(o) = c, /' (00) = 1, 7 



i/(x) = 7 log x + c, c = const. 
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